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Earth Currents Near a Monopole Antenna With 
Symmetrical Top Loading 

James R. Wait 

Expressions for the fields are developed for a vertical ground-based monopole with a cone or disk located 
at the top of the antenna to simulate umbrella top loading. The current distribution on the structure is 
assumed. Using spherical-wave functions, the magnetic-field distribution on the ground plane near the base 
of the antenna is computed and illustrated by graphs. For the case where the antenna is electrically small, 
the currents flowing on the cone or disk are shown to contribute only slightly to the total field. 

1. Introduction 

Transmitting antennas at low radiofrequencies usually take the form of a vertical current- 
carrying wire loaded by extensive overhead wires which sometimes take the shape of an um- 
brella. At the risk of oversimplification, a model for field calculations might be a vertical 
monopole with the umbrella being represented by a continuous cone. In other words, any 
asymmetry about the vertical axis is neglected. It is the purpose of the present note to inves- 
tigate the expected influence of the umbrella or cone loading on the fields produced by the 
antenna with an assumed current distribution. 

2. List of Principal Symbols 

r 1} #i, and <f>i = spherical coordinates of a variable point on the conical sheet. 
a = distance from the apex to the rim of the cone. 
J = current at the base of the vertical member. 
F( z) = current distribution along the vertical member. 

/&=height of the antenna or length of the vertical member. 
TT z =z component of the Hertz vector of the currents on the vertical member. 
K=the slope of the current distribution on the vertical member (assumed 
constant) . 
/0'i)=the current distribution along a generator of the cone. 
A Q = coefficient for the power series expansion of f(r { ). 
s' — a variable point on the vertical member (the integration variable). 
p, 4>, 2= cylindrical coordinates of the observation point P. 

H v t = tangential magnetic field on the ground plane resulting from the currents 
on the vertical member, 
lip and Il2=the p and z components of the Hertz vector of the currents on the cone. 
R= distance from the variable point Q on the cone to the observation 
point P. 
r, 6, z= spherical coordinates of the point P. 
P™(cos 0)=associated Legendre polynomial of argument cos 6 and orders n and m. 
jn(kr) = spherical Bessel function of the first type of argument kr and argument n. 
hn\kr) = spherical Hankel function of the second kind of argument kr and order n. 
T n {r u r)=jn(kr x )K%\kr). 

H d = tangential magnetic field on the ground plane resulting from the currents 
on the disk or the cone of angle ^ = 90°. 
A n (r ) } B n (r Q ), C n (r )= coefficients for the spherical-harmonic expansion of H d . 

r and 6 = spherical coordinates of a variable point on the ground plane for the 
origin at in figure 1, r =(p 2J rh 2 ) and o =&rctan (p/h). 
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j v =—aH v t /I . . . = earth current (dimensionless) on the ground plane resulting from the 

currents on the vertical member. 
j d =—aH d t /I . . . = earth current (dimensionless) on the ground plane resulting from the 
currents on the disk. 
r, 0=new spherical coordinates for describing the radiation pattern (see 
%• 3). 
Jo(S), Ji(S), and J 2 (S) are cylindrical Bessel functions of the first type of orders 1, 2, and 3 of 
the argument S. 

3. Some Preliminary Considerations 

The top-loaded antenna in its idealized form is shown in figure 1. It consists of a vertical 
current-carrying wire extending from the point 0, at z==0, along the axis of a cylindrical 
coordinate S3^stem to the ground surface at z—h. The point is also the apex of a cone 
which is generated by a line of length, a, at an angle X to the z axis. With respect to a spherical 
coordinate system with the origin at 0, the coordinates of the bottom rim of the cone are {a, 
6i, <j>) and those of a variable point Q on the cone are (r b U fa). The upward current on the 
central wire filament is taken to be I F(z) where I is the current at the apex, i.e., F(Q) = 1, 
and I F(h) is the current at the base. The outward current on an element dr x dfa on the 
conical surface is then (I /2T)dfaf(ri) where /(0) = 1 and/(a) = 0. Because of symmetry the 
current on the conical sheet has only a component in the direction of increasing r x . 




Cone Image 



Figure 1. Coordinate system for monopole with cone leading m 



248 



The current distributions, characterized by the functions F(z) and/(ri), in general, have 
a very complicated dependence on the antenna dimensions. In the absence of a top loading 
and where the vertical member can be approximated by a circular cylinder, current distribu- 
tions can be expressed in terms of an expansion parameter which is proportional to log (h/2irc) 
where c is the radius of the conductor [l]. 1 When h is less than one-half wavelength, the 
distribution is very close to being sinusoidal so long as h is greater than about 100 c. In fact, 
at very low frequencies where h becomes small compared with the wavelength, the current is 
then almost a linear function of height. It would not be expected that the cone would modify 
this conclusion. Therefore 

F(z)~I (l-Kz/h), (1) 

where K is a constant which depends on the extent of the top loading. For example, in the 
unloaded case, K=l, and in the ideally loaded case, K=0, corresponding to a constant current 
on the vertical member. In most applications where top loading is employed it is known 
on experimental grounds [3] that K~0. This is particularly so for high-power installations 
at low radiofrequencies where the effective dimension of a is comparable to or greater than 
h, and 0i is near 90°. 

Information regarding the current distribution function J(ri) for the conical sheet is 
apparently not available in the literature. However, if each individual wire making up the 
cone acts independently, it could be expected that the total current flowing towards the rim 
of the cone would vary in a linear manner with r, if a is small compared to the wavelength. 
On the other hand, when the cone degenerates into a disk (i.e., ^«90°) and, if h is somewhat 
smaller than a, it can be expected the total current flowing towards the rim will vary with 
the square of r x . This conclusion is based on the idea that the shunt displacement current 
density between the circular disk and the ground is a constant if fringing effects at the edges 
are neglected. A plausible assumed distribution would be 



■w-s®'^ ® 



where the summation could be over q=0 and 1 and A =l, Ai= — 1 for the linear current 
approximation; whereas g=0 and 2 with A = 1 and A 2 = — l for the squared current approxim.i - 
tion. It would be possible to employ more elaborate forms oij(r x ) by utilizing more than two 
values of q. A prime requirement is always 

S4«=0, (3) 

Q 

since/(a) = 0. 

4. Field of the Vertical Member 

-> 
The Hertz vector IT at the point P(p y 0, z) due to the current flowing only on the vertical 
-» 
member is given by IP=(0, 0, H v z ), where 

tt* 7 » / P exp [-ik[(z-zT+m m ,w, f exp [-ik[(2h-z'-zy+ P ^] V( ,, f] „, \ 

(4) 
when &=27r/wavelength and e =8.854X10~ 12 . The fields are then given by 

EI=(p+g)n?, and #;=-w> a 



#;=££' E'MV+^W, and #£= -ieoa, — *• (5) 



1 Figures in brackets indicate the literature references at the end of this paper. 
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The second integral in eq (4) can be identified with the image of the vertical member below 
the plane z=h. It can be verified that E v p = at z=h, which is the required boundary condi- 
tion for the ideally conducting ground. 

The tangential magnetic field, H t , on the ground plane is then explicitly given by 

rr, A A f F(Z') exp {-ik[(z'-hy + P*P} . , ,., 

n '- 2rbpJ ~ [( 3 '_A)2 +p2] M az - W 

Considering the important case where /3A<<1 and F(z')~l, it follows that 

g».-_. <-°>>I 9 fT 1 ■ ik - \ dz , 



[(?q^+^ tan_1 ®} 



(7) 



2xp |_(p 2 +#)> 
This can be expanded in a power series in (h/p), so that 

^--^(^(^^-(X^^---^^^^^® 5 ' ©' etc -} (8) 

The first term of this expansion corresponds to the field of a hertzian dipole radiator. The 
additional terms can be regarded corrections to account for the finite value of the antenna 
height and can be called multipole terms. Sufficiently close to the base, such that &p<<l, 
the following expression will suffice 

ttv ±o_ h , Q v 

which is a quasi-static approximation. A check on this result is obtained by noting that 

which could be predicted from Ampere's law. 

5. Field of the Conical Currents 

From an argument based on symmetry, it can be seen that the Hertz vector IT at P, due to 
the currents flowing in the cone, has a vanishing </> component. The p and z components, 
however, are, in general, nonvanishing. These can be expressed as an integration over the 
area of the cone, such that, 

nc = (/ /27r)sin^f 2 ^ f a ^ _ e~ ikR 



47ri€ co 
and 



I P2T f*a p -ikR 

Jo ^J/ riC0S ^~ir^ ri )> ( 10 ) 



where R 2 =r 2J r r\—2rr l cos 12 and cos fi=cos 6 cos ^ + sin 8 sin d 1 cos fa. Now the factor l/R is a 
solution of Laplace's equation and can be written in terms of spherical harmonics, as is well 
known [2], 

4=± ± < m feSi P « < C0S 6) P » ( C0S ^ C0S M ^ iKC! 5 ^ ^ 12) 

where e = 1, € m — 2 if (m^O), and where the P% are associated Legendre polynomials. The factor 
e~ ilcR /R, on the other hand, is a solution of the wave equation and an expansion for it in terms 
of spherical wave functions can be expected. Such an expansion does not seem to be available; 
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it can be obtained, somewhat heuristically, by noting that it should reduce to eq (12) for kr 
and kr x small, and should give rise to outgoing wave at infinity. Such an expansion is 

££--*± ± e m ^(2.+l)P»(cos *)P?(cos*0 ™n*I*™££\> r >« (13) 
n n =o m=0 {n+m)l {^jnikrjhn^ikr^ir^r, 

where jn(kr{) is a spherical Bessel function and h™ (kr) is a spherical Hankel function of the 
second kind. These are defined by 

>(*)= ( -i)»<^)W> ^ 



and 



The integration, with respect to </> h in eq (10) and (11) requires the evaluation of 

~*2ir — ikR 
LIT J o 1\ 

and 



(16) 



which can now be rewritten 



1 C 2lr e~ ikR 
I\ (rfi) =^: I — g- cos frdfa 

i r 27r p~ ikR 

'■m-sJ. -ir^> (17) 

ma—* § §rr] Pi(cos emcos e dkl:'r]\ T>r (18) 

and 

L(r,6) = -ik ± (2^+1)7^(008 0)P«(cos d 1 ) i l^'X'. '^ (19) 

where T n (r u r)=jn{kr x )h^ (kr). 

The circumferential field H^ due to the conical currents is given by 

ff$=ie « \-^ cos 0~^J> sin — gp- sin 0-^ cos 0j, (21) 

where 

ieoWn c = A^l f"/(n)/i(r,»)rfr lf (22) 

and 

teowIIS-^^^jVwii^dri. (23) 

To satisfy the boundary conditions at the ground plane, z=h, the fields of the mirror image of 
the cone must be added to H^. For these fields Ii(r,$) must be replaced by — I Y {r' ,6') and 
Mr,0) by +I 2 (r',d') where $' = M-d, r' = [(2h-z) 2 + P 2 f, and cos M=[(r)*+(r') 2 —4h*]/2rr'. 

The preceding equations, along with eq (7) and (8), constitute the complete formal solution 
for the fields of a cone-loaded monopole for a specified current distribution on the structure. 

6. Disk-Loaded Monopole 

When the angle 6 X approaches 90°, the cone degenerates into a disk. It is this case that 
will be considered in some detail here. Of particular interest is the behavior of the currents 
excited on the ground plane in the vicinity of the disk. In fact, for vlf antennas these near- 
zone currents have, associated with them, considerable ohmic losses due to the imperfect conduc- 
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tivity of the soil [4]. The surface current density is by definition numerically equal to the 
tangential magnetic field H d on the ground plane. For the case kr<^<Cl, it follows that 



and 



where 



H d =^ fj [4 n (r )Pi +1 (cos 6 ) +JB n (r )Pi_ 1 (cos d )]Pi(0) for r <a (24) 

H1=^ S CWPin (cos 0„)Pi(O) for r >a, (25) 

^ w= -r^HK^" n ' (26) 

^W=J ro ^r-^i> (27) 

^w— r*H^* ,; (28) 

and r = (p 2 +A 2 ) 1/2 and o =arc tan (p/A). 
Setting 

f(ri)=T,Mrila)«, (29) 

it is a simple matter to arrive at 

Cfr> y, 1 (ah)'*' A, 

Setting ^4 =1, A 2 = — 1 and the other ^4/s equal to zero, the squared current distribution, 
f{ri) = l — {r\la 2 ) ) is attained. Numerical values of the normalized current, j d = —aH?/I 
and j v =—aHt/I are shown plotted in figure 2 as a function of p/a for the squared current 
distribution on the disk and a constant current on the vertical member. The values of the 
ratio h/a are shown on the curves. The resultant current j on the ground plane which is 
radially directed is given simply by j d -\-j v in the near-zone region where p is small compared 
to the wavelength. 

It is of interest to examine the radiation field such that ?\>^> wavelength. The spherical 
Hankel functions can then be well approximated by the first term of their asymptotic expan- 
sion. The radiation field H d of the disk by itself is then given b}^ 

H " ~ - 1 ^£^ cos e S S4I P » (° os V P » (°> P M*ri)/(ri)<fri. (33) 

47TT n = in\Jl-\-L) Jrj=0 

The radiation of the disk and its image is obtained by adding Hi to the preceding equation 
which is of the same form as —H^ with ir— 6 replacing 6, and r—2hk cos replacing r. The 
summation in the preceding can be converted to a Bessel function by first starting with the 
addition theorem; 

cxp (ikn sin 6 cos 4>) = H S «- L lJu (2«+l)P?(cos 0)PJ(O) cos m<j> jnihrj , (34) 
m=o«=o (ri,-|-mj! 
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Distance From Base _ _P_ 
Radius of Disk "a 



Figure 2. Earth current distribution below a disk-loaded monopole. 
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and then multiplying* each side by cos </> and integrating with respect to <t> from to 2w. This 
leads to the interesting relation 



iJi(*ri sin9)=S §?Tn ^(oos (?) Pi(0) j»(*r.) . 



Therefore 



rj d kl e- ikr 
H%~— cos 



P Jin) Mkr, sin 9)dr, 

Jr,=0 



(35) 



(36) 



subject, of course, to kr» 1. The squared-current assumption, f(r 1 ) = l — (rf/a 2 ) is now intro- 
duced in the above which fortunately enables the integration to be carried out to yield 






cos 9 T(S), 



(37) 



where T(S) = [l—J (S)—J 2 (S)]S- 1 with S=ka sin 6. The field of disk plus its image can then 
be written 



/ « (1 P ^ — — — 

H$+gf« ° 9 - cos 6> i sin (M cos (9) T(fc sin 0), 



(38) 



in terms of the spherical coordinates r, d with the origin at the base of the vertical member as 
indicated in figure 3. For the usual case where ka sin 0<<1, it follows from the power series 
expansions of the Bessel functions J and J 2 , that 



1 psin 6) « ^ j 



(39) 




Figure 3. Coordinate system for describing radiation patte 
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where terms in (ka sin 0) 3 , etc. have been neglected. Noting that the radiation field for the 
vertical member in an /;, 6 coordinate system can be written 

Hl= %ICl f- sing, (40) 

it follows thai the resultant radiation field is given simply by 

g,~^(l- (fa)2 8 COS2 *) (41) 

neglecting terms of order (ka) 4 . 

7. Conclusion 

The results indicate that the earth currents near a monopole antenna are not significantly 
modified by the presence of a symmetrical top loading. Therefore, in any calculation [4] of 
the energy absorbed in the ground system, it is usually sufficient to consider only the contribu- 
tion from the currents in the vertical members (downleads). Furthermore, the radiation pattern 
of a loaded monopole at low radiofrequencies is essentially a dipole pattern being modified to 
small order by the currents in the top-loading structure. The effect of unsymmetrical top 
loading is considered in a sequel to this paper [5]. 



The author thanks Mrs. A. Murphy for carrying out the calculations and W. W. Brown, 
A. D. Watt, and H. A. Wheeler for informative discussions. 
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